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Abstract — Economic and Emission Dispatch (EED) 
represents one of the most important problem in Power system 
engineering. The major part of the power generation is due to 
fossil fired plants and their emission contribution cannot be 
neglected. This emission has to be reduced by minor change in 
the dispatch problem which results in increase in the fuel cost. 
There is a best compromise required in minimizing the emission 
and cost on pareto-optimum curve [6]. A direct method for cost 
optimization including losses was published by 2 authors in [8] 
& [7]. In this paper, direct method formulae taken from above 2 
journals & described for emission and combined fuel cost and 
emission optimization for achieving the best compromise. This 
method needs the calculation of combined coefficients obtained 
by applying the concept of proposed total deviation approach. 
This paper presents derivation & best solution for combined 
optimization of emission and cost with and without transmission 
losses for optimal solutions. The effectiveness of the proposed 
method is demonstrated by considering IEEE 30 bus system. 
The results for direct methods have been compared. The results 
confirm the potential and effectiveness of direct method 
compare to conventional lambda iteration technique and also 
verified using load flow. 


Index Terms — Combined Optimization, Cost & Emission 
co-efficient, Direct method, Economic and Emission Dispatch 
(EED). 


I. INTRODUCTION 

Generally the coal used in thermal generation is of poor 
quality and high ash. Sulphur dioxide (S0 2 ) and Oxides of 
Nitrogen (NO x ) are the major emissions from thermal plants 
due to the combustion of coal which will cause ill effects on 
human beings as well as animals [3].The NO x emission is 
required to be reduced by 2 million tons/year from 1980 level 
[4]. 

Minimizing operating cost can no longer be the only criterion 
for dispatching electric power due to increasing concern over 
the environmental consideration [5].The quadratic fuel cost 
function is well accepted for optimization problem. The direct 
approach has been presented for optimum scheduling in terms 
of total system demand P D [1]. 

Total fuel cost in $/hr 
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Fc = y F„(P„) 

L=1 

F n - Total fuel cost of n th generator in ($/hr) 


P n - generator power output of n th generator in (p.u) 

n 

Jl r = ^ 1 /2a„ 

: = l 
n 

Y = ^b n /2a n 



n- Corresponding generator (1,2,..N) 

N- no. of generator. 

a n , b n c n - fuel cost coefficients 

The closed form expression for calculation of 
optimum fuel cost has been given in terms of total system 
demand P D . All these closed form expression depends on the 
calculation of X, Y, Z calculated from the cost coefficients of 
the generation [2]. Optimum Generation Schedule 
P n * = ( P D +Y-b i X)/2a i (p.u) 

Optimum fuel cost 

Fc* = [(P D +Y) 2 /2 X]-z ($/hr). 


Co-ordination equation for determining the optimum 
generation schedule with transmission losses is 
dF„ _ dP L 


+ ^-7TT- = * 


dP 7 

Where, 


( 1 ) 




PR P 


Subjected to, 


x 


Pn ~P L ~Pr= 0 


dF n V 

-.(1 + 2 £iW*n) 


rr L = l 


& also Neglecting the off-diagonal elements of Bmn assume 
the form, 

A n Pn + B n P n + C n = X 


Where, 

4 = 40An. ^ = 2MVU C n = K 

An, B n , C n - functions of cost function coefficients and 
elements of 4 [Bmn] ’ , X ’ obtained from, 

d)? "I - + y\ "I" $ ~ Pr — 0 

Where a,p,y,d are expressed as follows, 
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Optimal generation scheduling is calculated by 
B n , [Bl+AA n {A-C n )f/2 


Pn ~ 2 A„ + 


2 A„ 


In [7], cubic should be solved. 

In [8], there are no such assumptions to found the direct 
formulae,Ignoring the off-diagonal elements of in 

equation (1) is equivalent to assuming „P m = 0" for „m#n \ This 
introduces errors in computing P^’ and hence in „ ’.Also 
these errors can be minimized by assuming „P m = P n which 
modifies the values of [8] 




- £ 


m=l 


For any values of ’ 


dF n 

— = z an p n +b ri 


dP L 


dP r 


p 

n n r n 


Substitute above equation in (1) &get, 

u (n). X 2 + v(n). X + w(n) 


i 

n=i 


Where 


0+Q(n)) 2 


Ll(?7) = 


-P* =0 


4 i? 


v(n) = 


a. 


2 B‘ 


2 a n b n + B, 

w{n) = - 


nn u n 


4Bi 


^nn 

Optimal generation scheduling is calculated by 

„ _ A-b„ 


^ C^Tl 




In [8], quadratic equation should be solved. 


In this paper, [8] is taken to solve emission and 
combined fuel cost and emission optimization, since all are 
quadratic functions. The u,v,w,Q constants depends on the 
coefficients of respective quadratic functions and Bnn. The 
minimization of emission and cost optimization problem 
requires combined coefficients. These coefficients are 
obtained by applying the concept of proposed total deviation 
approach. This approach is the minimization of the sum of 
cost deviation from the optimum cost and emission deviation 
from the optimum emission. Transmission losses are included 
by using Bmn coefficients. 


The effectiveness of the proposed direct approach in EED 
problem with & without transmission losses are demonstrated 


in IEEE 30 bus system. The optimal solution results are 
compared with Direct method technique [7] and also with 
conventional lambda iteration technique. 


II. PROBLEM STATEMENT 

Consider a system where 4 n’ be the total number of 
generators. The EED problem is to minimize the two 
competitive function total fuel cost and emission while 
satisfying several equality and inequality constraints [4]. The 
problem is formulated as follows: 


A. Objective Function 

Minimization of cost: For an 4 n ’ plant power system, the 
cost curves are expressed by a quadratic equation 

F(P n )=a n P n 2 +b n P n +c n $/hr ;n=l,2...N (2) 


Total fuel cost in $/hr 

fc = 2 

Equation (3) To be minimized 


(3) 


Direct expression for optimum fuel cost [2] 

Fc* = [(P D +Y) 2 /2X]-Z ($/hr) (4) 


Optimum scheduling [1] 


P n * = ( P D +Y-b n X)/2a n (MW);n=l,2,... N (5) 

Minimization of emission: Total emission of 
atmospheric pollutants such as Sulphur oxides (Sox) 
and Nitrogen oxides (NOx) caused by the fossil fired 
thermal generation expressed by a quadratic equation 
F(P fi )=Ae n P n 2 +Be n P fi +Ce n ton/hr ;n= 1,2.. .N (6) 

Total fuel cost in ton/hr 

Fe = E n (?) 

Equation (7) To be minimized 
Direct expression for optimum emission 
Fe *= [(P D +Y) 2 /2X]-Z (ton/hr) (8) 

Where X, Y, Z are calculated using A, B, C 
Optimum scheduling [1] 

P G i* = ( Pd +Y-Be n X)/2Ae n (MW); i=l, 2.. .n (9) 


B. Objective constraints 

Generation capacity constraint : Real power output 
of each generator with in lower and upper limits 
expressed as 

P n min < P n < P n max ;n=l, 2....N 
Power balance constraint: the sum of total system 
demand P D and transmission losses P L equals to total 
real power generation. Hence, 

P a =Y, P *~ P L 

M =1 

C. Problem Formulation(EED) 

Multi-objective problem is formulated into single 
objective problem as 
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Minimize [Fc, Fe] 

Subject to: P n min < P n < P n max ;n=l, 2.. ..N 

n 


= 


Y 

n-n = "I 




III. Principle of combined Optimization 

Multi objective problem is formulated into single objective 
problem in terms of percentage deviation, obtained from the 
individual optimum solution Fc , Fe . 

• Cost deviation function obtained from the individual 

optimum cost in percentage : 

Dc = ^^-**100 = (— — l)* 100 (10) 

Fc* \F c* / v 7 

• Emission deviation function obtained from the 

individual optimum emission in percentage: 

De =LlzIll * ioo = (—- {]*■ 100 (11) 

Total deviation function: 

Td = Dc+De (12) 

From the equations (1), (5), (9) & (10) in (11). 


ll-J{Pn)~ F(?n) • Z.J(Pn)-E(Pn)’ 


F(Pn)‘ 


E(Pn)’ 


I^F(Pn) ] (MPn) 


F(Pn)- 


-1 


£(Pn)* 


i; = 1 F(Pn)l (g' sl E(Pn) 




F(Pn) ■ j ( f(Pn)» 

"{ F(Pn)‘ )*( E(Pn)‘ J 

-f [ a A h A _L_|+ yf k A B A _i_L 

L 1 F(Pn) * F(Pn) * P(Pm)‘] Z- [ £(Prt) * f(Pn)- F(Pn) • J ' 

n=l N ' .*i=l N 

“Xl*" L=*(Pn) * * E(Pn) * ] + Pk IfCPm) - ~ E(Pn ) * ] + L^Pri) * * E(Pn) * ] “ 2 } 

n=l 

= a n P n +Pn P n+Yn (13) 


Total deviation function coefficients are, 


n Fc* Fe* 

(14) 

rs _ bfl _ 

Fc* Fe* 

(15) 

r n Fee Fes 

(16) 


The transmission losses are included by ‘Bmn ’ coefficients. 

IV. PROPOSED DIRECT METHOD 

The closed form expressions depends on the total system 
demand P D are used in all the three cases of optimization 
(i)Cost (ii)Emission and (iii)Combined cost and emission. 
The approach is described in a common flowchart fig.l. 

i. Suboptimal solution: 

Co-ordination equation is taken as 


dF n 

dP u 


= A 


In this equation incremental transmission loss is changed at 
the equal incremental production cost value. 


V. RESULTS 

A. Direct method results 

In order to validate the proposed method, the emission and 
economic dispatch problem was solved for IEEE 30 bus 
system shown in fig.3 and results are presented in this 
section.There are 6-Generators.The total system demand is 
283.4MW. The values of fuel cost coefficients and emission 
coefficients are given in Table I. MATLAB program was 
developed to perform economic and emission dispatch (EED) 
problem. 

Results are obtained for Cost, Emission and Combined fuel 
cost and emission optimization for both with and without 
losses for optimal and suboptimal solutions. Results obtained 
demonstrate its effectiveness in finding optimal solution and 
compared with conventional lambda iteration technique. The 
computational result of the proposed method for the EED 
problem, where as total fuel cost is minimized as shown in 
Table II. From Table.II it is clear that fuel cost is minimized in 
EED problem and the corresponding emission is 
higher.Results of best emission are shown in Table III. Here 
the emission is minimized and the total fuel cost is higher. 
The fuel cost and emission are compromised in the 
optimization of total deviation function. Results of best 
combined fuel cost and emission is shown in Table IV. 

Graph 1. Clearly shows that, when fuel cost is minimized and 
emission is higher. When emission is minimized and fuel cost 
production is higher. In the combined optimization both fuel 
cost and emission are compromised. It gives a best solution in 
the pareto-optimum curve. 

From the best solution, if the fuel cost deviation is reduced, 
the deviation of emission will increase and vice versa. This 
optimal solution is rechecked by substituting generator power 
in load flow and from the result,it is found that same 
transmission loss is obtained. 

B. Comparsion of Optimal solutions 

Case 1: EED by lambda iteration technique 

The total generation costs and emission of easel are shown 
in Table V. The results show that the lambda iteration 
technique consumes least computing time when compared to 
other methods in terms no. of iteration is high and time taken 
also high. The total fuel cost calculated by lambda iteration 
method for demand of 2.834p.u is 616.64$/hr. 

Case 2: EED by Direct method [7] 

The total generation costs and emission of case 2 are shown 
in Table V. in this method there several approximations made 
in co-ordination equation and non-diagonal elements get 
eliminated. So the total fuel cost gets deviated. The total fuel 
cost calculated by direct method [7] method for demand of 
2.834p.u is 628.819$/hr. 

Case 3: EED by Direct method [8] 

The total generation costs and emission of case 3 are shown 
in Table V. in this method there is no such approximation 
made. Non diagonal terms also included. The total fuel cost 
obtained by this method is very nearly equal to lambda 
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iterative technique. The total fuel cost calculated by direct 
method [8] method for demand of 2.834p.u is 615.186 $/hr 
This method results is verified by the optimum generator 
power (P G ) is substitute into load flow get the same loss which 
is equal to loss obtained by this proposed method. 

The results obtained by the three methods for a 30 bus 
system are compared and it is shown as Graph 2. It shows the 
merit of direct method. The fuel cost of direct method [7] is 
equal to lambda iteration technique. Computation time is also 
very less in the direct method. The Graph 3 shows the 
comparison of optimal and suboptimal solutions for 
combined fuel cost and emission optimization. It clearly 
shows the merits of optimal solutions. The sub optimal fuel 
cost is nearly equal to optimal solutions. The percentage 
deviation between optimal & sub optimal is only 0.06899%. 


[4.] M.A Abido ‘Environmental/Economic Power Dispatch Using Multi 
objective Evolutionary Algorithms’ IEEE Transaction on power 
system,Vol.l8 No.4, November 2003, p.1529. 

[5.] Maheswarapu sydulu(1999),’A very fast and effective non-iterative X- 
logic based algorithm for economic dispatch of thermal units’, 
Proceedings of EEE Tencon. 
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generation schedule”Elsevier-journal vol-ll,pp.147-152,july 1986. 
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VI. Conclusion 

This paper presents the direct approach for combined 
emission and cost optimization for thermal power plants. The 
proposed method derived from basic work presented for 
obtaining closed form expression (in term of P D ) [8]. In order 
to convert the multi-objective (EED problem) into single 
objective problem, Total deviation approach is applied. The 
total deviation is the sum of cost deviation and emission 
deviation. This approach helps for arriving the best balanced 
compromised optimization solution for any multi-objective 
problem. Here optimal and suboptimal solutions are 
compared. Transmission losses are included. In optimal 
solutions losses are also included by direct approach. In order 
to include the transmission losses in sub optimal solutions, 
Two-Fold method is applied for feasible operating point of 
the system in EED. Combined optimization gives best 
solution by the proposed approach. The following 
conculsions are drawn from the proposed work, 

1. Shooted for all Real time application. 

2. Optimal and suboptimal solutions are not having 

much difference (ie., 5% when network loss is less 

than 25%). 

3. No iterative technique. 

4. Using this method Multi objective is converted into 

single objective problem. 

5. Easiest technique to find the best feasible global 

optimum. 

6. Transmission losses are included. 

7. The solution of combined cost and emission 

optimization is the best compromised. The solution 

is arrived by using the closed form expressions. 

8. Computational time is very less, since the closed form 

expressions are used for obtaining the solution. 
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